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Abstract 

We exhibit examples of sub-Riemannian metrics with integrable geodesic 
flows and positive topological entropy. 

Introduction 

Consider a distribution on a manifold Af™, i.e. subbundle of the tangent 
bundle 11 C TM. Non-holonomic Riemannian metric is a Riemannian metric 
g S S^IV on this bundle. We call the pair (11,(7) suh- Riemannian structure. A 
curve 7 : [0, 1] M is called horizontal if 7 is a section of 11. We denote the 
space of horizontal curves joining x to y hy 1-L{x, y). A theorem of Rashevsky- 
Chow states that if M is connected and 11 is completely non-holonomic 

then 'H{x, y) is always non-empty. By completely non-holonomic we mean distri- 
bution n, such that the module 'D'^'^ of order < N self-commutators (of various 
kinds) of sections of 11 is equal to the module V{M) of all vector fields for some 
big TV. From now on we consider only completely non-holonomic distributions. 

For horizontal curves we calculate its length lg{"f) = Jq \\j\\gdt and this 
produces sub-Riemannian distance (metric) on M by 

dg{x,y)^ inf lg{-f). 

-feH(x,y) 

A curve 7 € 7i is called geodesic if it realizes the minimum sub-distance for any 
two of its close points. The description of the most geodesies (normal ones) 
is given by the Euler-Lagrange variational principle. There is a Ham iltonian 



reformulation of this principle, due to Pontrjagin and co-authors PBGM |, which 
allows to consider the geodesic flow as the usual Hamiltonian flow on T*M . 
There appear occasionally geodesies of different kind - abnormals - which are 
not governed by the Pontrjagin principle for 7, but depend on the distribution 11 
only. However if we consider contact distributions H, i.e. distributions such that 
for any non-zero section a of the bundle Ann(H) C T*M we have a A (da)" 7^ 
for m = 2n-\-l (in particular m = dim M is odd) , then all geodesies are normal. 

As in the standard theory of geodesies we say the metric g is integrable if 
the Hamiltonian flow of this metric is integrable on T*M in the Liouville sense. 
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i.e. there are a.e. functionally independent integrals Ii = H, I2, ■ ■ ■ , In which 



Poisson-commute {7^,/;} = (see [BF] for a review of methods and problems). 



For any diffeomorphism / : Af — > Af the suspension MxfS^ is the quotient 
of M X R by the equivalence relation {x,t) ^ (f{x),t + 1). In particular for 
any matrix A G SL„(Z) we have a linear automorphism A : T" T" of the 



torus and hence the manifold Ma — T"xaS^. In the paper |BT1| a series 
of examples of (Riemannian) geodesic flows on Ma was constructed with the 
property that they are C°°-integrable, C"^- (analytically) nonintegrable and have 
the topological entropy htop > 0. 

We show that similar effect takes place in sub-Riemannian geometry too. 

Theorem 1 . For every A e SL2 (Z) there are a contact distribution H on 
Ma = T^XaS^ and a non-holonomic Riemannian structure g on it such that 
the sub-Riemannian geodesic flow is C°° -integrable. 

In the case of real eigenvalues for A ^ ±E the flow is not integrable with 
geometrically simple set of integrals. Moreover the topological entropy is posi- 
tive in semi-simple non-trivial case, htop = max {In |Ai |, In IA2I}, and vanish 

otherwise, htop — 0. 

The geometrically simple set of integrals includes analytic polynomial by mo- 
menta integrals and is defined in §^ (r emark that we consider only integrability 



with Liouville tori, not cylinders, see §1 (A) ) 



Note that integrability of sub-Riemannian flow is a more complicated (and 
fascinating) fact since the description of the geodesies is completely nontrivial 



even locally (see [ VG| for a picture of sub-Riemannian wave front). 

In particular locally over a point x of the manifold Riemannian geodesic flow 
in t'^/{Ux) is integrable, while in non-holonomic case it is no longer true. On 
the other hand integrability does not imply good behaviour in other respects: 
The well-known Martinet case (see ex.^) is C"^-integrable, but small geodesic 



balls are not sub-analytic ( |AS|] ) and there are abnormal geodesies. 

Remark 1 . One can easily generalize the above examples to the examples of 
integrable sub-Riemannian geodesic flows on Ma for arbitrary A £ SL„(Z) with 
n > 2. There appear distributions of various ranks and kinds (they correspond 
to higher Vandermonde determinants instead 0/ 2 x 2 matrices given below). 
Though completely non-holonomic, they have very degenerate curvature, i.e. 
the map Qu : A^II v ^ TM/Il induced by commutators (in particular, 
the growth vector can be rather long). They never produce higher- dimensional 
contact examples. 

Every sub-Riemannian metric g on a manifold of dim = 2n + 1 produces 
canonically a transversal to the contact distribution 11 vector field - Reeb field 
Va (also denoted Ra). To see it we choose a e Ann 11 such that || daljj || = 1, 
where the norm is induced by the metric. In other words ((ia)"|]-[ = volg. This 
gives a uniquely for odd n and up to ±1 for even if the distribution 11 is oriented. 
The Reeb field is defined now uniquely by a(i^Q) = 1, i^^da = 0. 
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Theorem 2 . In addition within the construction of theorem ^if Tl is coori- 
ented (Sp A is contained in either R-)_ or C C) the topological entropy of the 
Reeb flow is htop{va) = 0. 

We introduce Riemannian metric g on M by the requirements gjjj = g, 

Va -Lg n, \\iya\\g = 1- 

Theorem 3. The geodesic flow of Riemannian metric g on Ma is C°°- 
integrahle, -nonintegrahle for real eigenvalues of A =/= ±E and has positive 
topological entropy for semi-simple nontrivial case. 

In fact the metric g coincides with the metric on Ma introduced by Bolsinov 
and Taimanov. The first example of C°°-integrable geodesic flow violating ob- 
structions for C"^-integr ability was found by Butler & and then it was included 



into a much bigger family of examples-suspensions in |BT1]. Their examples 
obviously generalize to Lorentian and other semi-geometries. We show the case 
of sub-Riemannian geometry (actually of any dim) can be considered as well. 



Remark 2. Note that examples from |BT1] are easily generalized to a larger 



group GL„(Z) = SL„(Z) x Z2. However the sub-Riemannian case is differ- 
ent: the manifold Ma with det A = — 1 does not admit any contact structure. 
Actually a contact structure on a 3- dimensional manifold gives a canonical ori- 
entation. 

We discuss some relations between g, g and h'a at the end of the paper. In 
§ 1^ we exibit Poisson action of such that entropy of every non-zero vector is 
non-zero. 

Let us note that theorem |l| provides the first example of C"-polynomial in 



momenta nonintegrable sub-Riemannian geodesic flow. In paper |MSS| it was 
presented only algebraically non-integrable sub-Riemannian flow on a Lie group 
of dimM = 6. In dimension 3 all left-invariant sub-geodesic flows on Lie groups 
are Liouville integrable (as we show in §^). 

Acknowledgements. I would like to thank J. P. Gauthier, A. V. Bolsinov, V. S. Matveev 
and O. Kozlovski for helpful discussions. I am indebted to A. B. Katok for the reference in 
My interest to sub-Riemannian geometry has grown up during a Singularity Theory workshop 
at Newton Institute (Cambridge), 2000. I would like to thank the Institute and especially 
V. V. Goryunov and V. M. Zakalyukin for the organization of this activity and hospitality. 



1 Sub-Riemannian geometry and Hamiltonian 

systems 

(A) Integrable Hamiltonian flows. 

Let (ly^", n) be a symplectic manifold. Hamiltonian vector field is the field 
r2-dual to an exact 1-form dH: il{XH, •) = dH. We denote the field Xh also by 
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sgrad-ff, since the above is similar to the definition of the usual gradient (and 
in order not to confuse with Ma above). 

Poisson brackets can be defined as {F, G} — VI{Xf,Xc)- The Hamiltonian 
system is called complete (or Liouville) integrable if there are additionally to 
Ii = H involutive integrals I2, ■ ■ ■ , In, {Ij, h} = 0, which are functionally in- 
dependent a.e. By Liouville theorem (Q) a full measure set W C W is then 
foliated by cylinders (tori in compact case), and each cylinder has a neighbor- 
hood with coordinates tp € T^-^ x and new / G R" s.t. ^ ^ dl A dip, 
{Ij = Cj}i<j<„ ~ j^n-r ^ jg^r ^-^^ ^y^q flow \s ip ~ '^{I)- Thc numbcr r (noncom- 
pactness rank) vanishes in many important examples. We will always suppose 
r = (even for noncompact W) and call this case Liouville integrahility . 

Let (M, g) be a Riemannian manifold. Then the geodesic flow ipt ■ TM 

r 

TM can be considered as Hamiltonian on T*AI ~ TM with the standard sym- 
plectic structure Q if we choose the Hamiltonian H ~ ^\\p\\^, p e T*M. We say 
g is integrable if the flow (ft is Liouville-integrable. 



(B) Pontrjagin maxiinum principle. 

Consider now non-holonomic case. We start with an arbitrary completely 
non-holonomic distribution H C TM. 

Let {T*M, 51) be the cotangent bundle equipped with the standard symplec- 
tic structure. The non-holonomic metric defines isomorphism jj^ : H* ^ H. 
Consider the inclusion i : H ^ TAI . Then we have vector bundles morphism 

defined as the following composition: 

T*M ^ H* ^ H A TM. 



Contrary to Riemannian situation this is not isomorphism. We have: Ker(^'g) = 
AnnH and CoKer(^'g) = v = TM/U. 

Define the Hamiltonian function H on T*M as the composition 

T*M ^ H* ^ H M. 



This function can be locally described as follows. Let ^1, . . . be some or- 
thonormal basis of vector fields tangent to H. Every vector field is a fiber-linear 
function on T*M. So we have H = ^ EtiL__ 

Thc Pontrjagin maximum principle ( ||PBGM |) states that trajectories of this 



vector fields in region {H > 0} projected to M are optimal for the corresponding 
variational problem. They are called (normal) geodesies. 

Example 1. Consider M = T^ — R'^/27rZ'^ with cyclic coordinates {(pi,ip2,^3)- 
Let a — siiiipidip2 + cosipidip^ be a contact form and H — Kera the cor- 
responding distribution. Let the metric be induced from the standard metric 
ds^ = '^d(pf on the torus. Then H — ^[p\ ~\- {cofnpiP2 — snvfip-if'] and the 
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Hamilton equations ip = dH / dp, p = —dH/dip have the form: 

2, N sin(2i^i) . sin(2(y9i) . s 

Vi=Pi, V2= cos {(pi)p2 ^ P3, = ^ P2 + sm {fl)pi, 

sin(2(j3i) ,99, , , 

Pi = ^ ^(H - Vi) + COs(2^i)p2P3, P2 =0, P3 = 0. 

T/ie?/ can be easily integrated. The Junctions I2 — P2, I3 — Ps are obviously 
integrals. This can be also checked via Poisson brackets: the functions Ii = H , 
I2, I3 o,re involutive: {Ij,Ik} = 0. The only singularities are Q and QD (so- 
called atoms A and C2). 



Example 2 . The geodesies on the Heisenberg group are given by the 2H = 
(pi + X3P2)^ +^3- Solving the Hamiltonian equation we see that geodesies are 
spirals in the direction of X2-axis, projecting to arbitrary (including radius 00) 
circles on the plane ]R^(a;i, 0:3). 

(C) Note on abnormals. 

In some situations there are optimal lines, which depend on the distribution 
n only. They are called abnormal geodesies and in generic case are described 
as follows. 

Consider the submanifold S ~ Ann 11 C T*M. It is given by equation H = 0. 
So contrary to the Riemannian case k — codim5 = dim 11 < m = dimAf. Let 
Om C T*M be the zero section and Sq = S \ Om be the complement. This So 
is a punctured cone over Om- 

Let k E 2Z. Then generically there is a hypersurface l]2n-fc-i ^ (maybe 
empty), where r2"~'^/^|^ = 0. The integral lines of Kcr(il|j;) projects to the 
abnormals. 

Let k E {2Z + 1}. Then generically kernels of ilj^^ are one-dimensional, so 
they integrates to lines projecting to the abnormals. 



This geometric description comes from the variational approach of [ Hsu 



Example 3 . Consider the Martinet case: H — Ker a for a ~ dy — z^dx on 
M?(x, y, z). Let the metric be lifted from the Euclidean plane M^(.t, z) (recall the 
abnormals do not depend on a choice of the metric). So 2H = {px + z'^py)'^ -{-p^. 
There is no flow on S = {H = 0} (as usual because H is quadratic). However 
there are geodesies. 

Actually, S — {px — —z'^py,pz = 0}. So the coordinates are {x,y,z,py). 
Since = dx A dpx + dy A dpy + dz A dpz we get il^ | ^ = 2zpydx A dy A dz A dpy . 
Since on So py ^ we conclude that = {z = 0} C with coordinates 
{x,y,Py). Now rijj; = dyAdpy, w/ienceKerfi = (dx) and the abnormal geodesies 
are given by {y ~ const, z = 0}. 
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Example 4. Consider the Engel distribution: Ann 11 = {dy — zdx,dz — wdx) , 
i.e. n — {dx + zdy + wdz, dw) on M^{x,y,z,w). The manifold S is given by 
{px = —zpy — wpz,Pw = 0} and has coordinates {x,y, z,w,py,pz). So = 

is equivalent to pz = 0. This gives the submanifold and rt\.^ = —pydx Adz ~ 
zdx A dpy + dy A dpy, whence Ker(f2|j^) — (dw)- This is so-called characteristic 
direction for the Engel distribution and the abnormals are just the corresponding 
integral curves. 

It follows from the above description that non-holonomic metrics on contact 
distributions have no abnormals. 



(D) Local integrability. 

Though in many respects sub-Riemannian structures are similar to Rieman- 
nian ones (e.g. ]ACG| ), as the previous subsection shows there are differences. 



Another discriminating aspect is Liouville integrability. 

Let us consider Euclidean metric locally in a neighborhood (x) . Symplec- 
tic reduction of the corresponding Hamiltonian flow on {H = 1/2} is T^^S"^^. 
By convexity theorem for arbitrary Riemannian metric and small e the sym- 
plectic reduction preserves differentiable type of T*S^~^. Moreover by Moser 
isotopy method applied to (open) neighborhood of the zero section we get that 
germs of reductions are symplectomorphic in Riemannian and linearized (Eu- 
clidean) cases. So one can find involute set of functions and integrate (by Liou- 
ville) the Riemannian geodesic flow on a neighborhood Ux C B^{x). 

In contrast even the differentiable type of symplectic reduction for sub- 
Riemannian structure in the simplest Heisenberg case can be changed under 
arbitrary small perturbation. That is because the geodesies are curved and 
being perturbed may leave/return the domain. 



2 Entropy of a dynamical system 

(A) DeGnition of entropies. 

Consider at first discrete time dynamical systems. Let M be a compact 
topological space with probability measure fj, and / : M ^ M be a homeomor- 
phism. Consider some partition ^ (up to zero measure) M — U^^ of our space 
by positive-measure sets with finite or countable number of indices. We form 
the new partition ~ ViLo / '"^j where ^ V 77 is a partition formed by the 

sets n 77/3. 



We define (|KH]) entropy of the partition by H{^) = — ^ 1^- and 
entropy of a preserving measure map / w.r.t. ^ by hu{f;^) = lim ^ — —. 

n — >oo n 

Then we define measure entropy of / by 

= sup hi,{f;^). 
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To define topological entropy we should change partition ^ by an open cover it 
of M = ^Uj, function H{^) by the function -/V(il) that is cardinality of minimal 
subcovering and then take sup over all open coverings. This quantity is subject 
to the following variational principle: 

hopif) = sup/ip(/), 

where the supremum is taken over all /-invariant Borcl probability measures on 
M. 

Another useful definition of the topological entropy goes for metric spaces 
{M,d). We define — max {f^Yd and denote by S{f,e,n) the minimal 

0<i<n— 1 

number of e-balls in dJ^-metric to cover M. Then 

/ztop(/) = lim IS ^J^^illll^ ^ Mm hm i^i^iZlflZ^. (1) 



This limit depends not on d but on the d-topology only (|KH]). 

Example 5. Let A e GL„(Z). Then we can define map A : T"- ~* T"^ , 
where T" = R"/Z". Let the measure be dfi — dxi A ... A dxn in the standard 
coordinates. Then 

htoj,{A) = hM) = 1^1^^ I' (2) 

Aa6Sp+(A) 

where Sp^(^) is the part of the spectrum outside the unit disk Di C C. 

If V € T>{M) is a vector field on a manifold M we define hto-p{v) = /itop(</'i), 
where tp* : il/ — s- M is the flow of v. We define the measure entropy similarly 
and note that for continuous time dynamical systems the entropies are also 
connected by the variational principle. 



(B) Compactness restriction. 

If M is noncompact, the previous definitions basically do not work. Still it 
is possible to define the entropy. 

The easiest case is as follows. Suppose there exists a compact exhaustion 
M ~ Uj^iKj with Kj C Kj+i and each Kj being compact and /-invariant. 
Then we define topological entropy by 

htopif) = lim /itop(/lK.) e M>o U {+oo}. 

Since the sequence /itop(/lx ) non-decreasing the limit exists and since it 
exists for every exhaustion {Kj}^ it does not depend on its choice. We call 
such case non- compactness of the first type. 
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As an example we consider Liouville integrable Hamiltonian system sgrad H 
with commuting integrals Ii = H, I2, ■ ■ ■ , In on M^" , for which generic intersec- 
tion {Ij = const j} is a union of tori. Then we obviously have invariant compact 
exhaustion and define /itop (sgrad i/). 

In the general situation we consider arbitrary compact exhaustion M — UKj 
(not necessarily /-invariant). Let us study the restriction of the metric to 
Kj. Since Kj is compact there is a finite covering by e-balls. Calculating the 
asymptotic of the minimal cardinality of such a cover we get similarly to (|l|) the 
quantity htop{f',Kj)- This sequence is non-decreasing by j and we set 

htopif) = hm htop{f;Kj). 

Again the limit does not depend on a choice of the compact exhaustion. But 
it can be infinite even for sufficiently smooth /. Let's call the situation non- 
compactness of the second type. 

Example 6. Consider a linear automorphism A : R" ~* R". Then the entropy 
htopiA) is given by the same formula as in example |^. 

Note that this example shows /itop(/~^) 7^ ^top(/) generally in non-compact 
case. However all properties of the entropy holds for non-compactness of the 
first type. 

The measure analogs of the above definitions are straightforward. 

(C) Lyapunov exponents. 

Let II • II be any norm on TM and / : Af ^ M be of class C^. Then we 
— ln||/^w|| 

can consider the map v i— > lim — . This map takes values Xi{^) 5: ■ ■ ■ 5: 

fc— >oo k 

Xn{x) at (almost) any point x € M", which are called Lyapunov exponents 
(multiple values repeat). They are defined similarly for the flows and do not 
depend on a choice of the norm defining the topology. 

Ruelle inequality says that if fi is /-invariant Borel probability measure, then 



M 



In particular if the norm can be chosen so that ||(<pt)*'u|| is constant (Lebesgue) 
a.e., then the entropy vanishes w.r.t. any Liouville (i.e. absolutely continuous) 
invariant probability measure. Therefore we get 

Example 7 . If a Hamiltonian system with Hamiltonian H is completely inte- 
grable, then for every Liouville invariant measure /i^(sgrad_ff) = 0. 

Moreover Pesin theorem states that the inequality above becomes equality 
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for C'^+'^-difFeomorphisms/flows and any Liouville measure. 

Since for integrable Hamiltonian system every measure p{Ii, ■ ■ ■ ,/n)|r2"| is 
invariant, where fl is the symplectic form and Ij are integrals, there are finite 
Liouville measures in this case (even if Af^" is non-compact). 



(D) Geodesic Bows. 

Let (M, g) be a compact Riemannian manifold. Let ipt : TM — > TM 
be the geodesic flow. Note that restriction of this flow to invariant compact 
submanifold TcM — {v G TM, \\v\\g — c} is conjugated to ipt\rp_^]yj. Hence we 
define entropy h{g) of the metric g to be the entropy of the last flow. 

Note that since the geodesic flow is Hamiltonian the (Liouville) measure 
entropy vanishes hf^{g) = for an integrable metric g due to example ^. The 
behavior of topological entropy as we will see in §^ can be more complicated. 

Example 8. Geodesic flow on the group S0(3) with left-invariant Riemannian 
metric has htop{g) = 0. 

Remark that canonical identification T*M ~ TM allows to consider the 
lifted flow (fit : T*M T*M. Its restriction to the hypersurface T^M = {H = 
1/2} determines the entropy, where H is the Hamiltonian of the geodesic flow 
(for non-compact M we should use the above modification of /itop)- 

Consider now sub-Ricmannian structure {M,Il,g). We again can consider 
geodesic flow and restrict it to the isoenergetic surface Q = {H = 1/2} C T*M, 
ifit '■ Q Q- This Q however is always non-compact, so that we should use the 
non-compact version of the entropy: htop(Ji,g) '■= htop^ftlq)- 

Example 9. For the group S0(3) the topological entropy of sub- Riemannian 
geodesic flow (see §y[) is zero. This follows from the study of bifurcational dia- 



gram ( also directly from the description of the geodesies in , VG / ) 



Note that if the Hamiltonian flow sgradi? is Liouville-integrable we have 
non-compactness of the first type, Lyapunov exponents vanish a.e. due to Liou- 
ville theorem and the (Liouville) measure entropy is zero. 



3 Commuting dynamical systems 

Let f,g:M—fMhe two commuting dynamical systems. They define an 
action of on M. If we have more commuting systems or commuting flows, 
then higher dimensional groups Z'= or K*^ act. There is a definition of entropy of 
such an action (§). However if it is non-trivial, then entropies of all dynamical 
systems given by 1-dimensional subgroups are infinite. We consider instead the 
case when the higher entropy is zero and the entropies of the generators of the 
group are finite. 
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Theorem 4 ([Hu|). Consider a 'L'^-action with generators f,g of class C^~^^ . 
Let h be either h^, where /i is any "Z/^ -invariant Borel probability measure on 
M, or htop and then we additionally suppose fg £ C°° . Then: 

h{fg)<h{f) + h{g). (3) 

Proof. Let us indicate the proof for a Liouville measure. Let xl i^) and 
Xj{x) be Lyapunov exponents for / and g respectively. Take xli^) < A < 
Xi+ii.^)- Then the map Fx :— A^^/, : TM — > TM commutes with and so 
{F\Y o g* — g* ° (FxY . Tending r — > ±00 we get y-invariance of the repeUing- 
expanding directions decomposition TxM = H~ © ■ 

Applying this to all A different from Lyapunov exponents and interchanging 
/ and g we deduce for x from a full measure subset in M the common decompo- 
sition TxM = (BHs{x) which is: measurable, invariant and such that on every 
subspace the Lyapunov exponents of both / and g are constant. There- 
fore the Lyapunov exponents oi f o g are sums of the corresponding Lyapunov 
exponents for / and g. 

The claim now follows from the Pcsin formula |(C), The equality in (||) 



is achieved if and only if we sum positive exponents for / and g in the above 
decomposition (SHi almost everywhere w.r.t. fi. □ 
Th e case of general Borel measure as well as the inequality for htop is done 



in |Hu| using Lyapunov charts (theorems B,C). 



Remark 3 . Without commutativity assumption the statement is false. Actually 
the Anosov automorphism A : T^ ^ T^ with ^4 = ( j can be decomposed 

A = BE* with B = But htop{A) = In ^ , while htop{B) = 

htop{B') = 0. 

Also the requirement that the dynamical systems are smooth is crucial: There 



are continuous -actions ([Pa]) with vanishing directional entropies for all ir- 



rational directions and non-zero for all rational. 



Corollary 1. Let I?' -dynamical system given by smooth commuting f,g be 
either defined on a compact manifold or have non- compactness of the first type. 
Let h denote h^ or htop- If h{g) ~ 0, then h{f g) = h{f). 

Proof. Actually h{fg) < h{f) and using h{g^^) ~ h{g) = we get the 
reverse inequality. □ 

Similar to (|^) formula holds for commuting flows ip*^ , (p!^ generated by vector 
fields ^,77: 

h{ipy^)<hi4) + hi^:^), (4) 

where h = /itop or h^ as in the theorem. 
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Therefore if R'^ acts on M with generators (p\^ , • ■ • , we define 

p{v) ^ h(^\' ■ ■ ■ ^l') for v = {h,...,tk)eR''. (5) 

The inequaUty (Q) can be seen now as triangle inequality and hence p is a 
pseudonorm on R*^ (provided non-compactness is no greater than the first type). 
This p however should not be continuous and can be degenerate (not norm). 

4 Integrability of sub-Riemannian metrics on 3- 
dimensional Lie groups 

Consider 3-dimensional Lie groups G. Every left-invariant structure is pre- 
sented on the Lie algebra Q. So we write the classification of left-invariant 
non-holonomic metrics on contact structures along Bianchi classification of 3- 



dimensional Lie algebras (| VG|). We identify proportional metrics (the geodesic 
flows are reparameterized). We present distribution H <Z S and the metric by 
an orthonormal basis ^i, 

1. Heisenberg algebra t)(3): [61,62] = 63, [6^,63] — and — ^k, k — 1,2. 

2. Solvable algebras: [61,62] = 0, 

(a) [61,63] = A161, [62,63] = A262 (Ai 7^ A2) and ^1 = 61 +62,^2 = 63. 

(b) [(61 + 162), 63)] = 6^*"^ (61 + 262) (0 < </? < tt) and ^1 = ei, $2 = 63. 

(c) [61,63] = 61 + 62, [62,63] = 62 and ^1 = 61,^2 = 63. 

3. Orthogonal algebra so(3): [61,62] = 63, [62, 63] = 61, [63, 61] = 62 and 
Ci = ei,^2 = cre2. 

4. Special linear algebra sl2(R): [61,62] — 63, [61, 63] = 26i,[62,63] = —262 
and two possibilities 

(a) Ci = 61,^2 = CT62; 

(b) 6 = 61 + 62,6 = crea- 

In three last cases the parameter cr G R+. 



One can integrate the equations of geodesies (see [ VG [ for a description in terms 



of semi-direct product). We consider Liouville integrability. 

Theorem 5. Non-holonomic geodesic flows of the above metrics on 3-dimensi- 
onal Lie groups are -integrahle. 

Proof. First note that any vector v € Q generates left- and right-invariant 
vector fields and i?„ on G, which can be considered as functions on T*G. 
We also denote by Ly,Ry e C°°{T*G) the invariant extensions of any function 
V e C°°{G*), not necessarily linear. Any two L„ and Poisson commute and 
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the other Poisson brackets are {Lt,,Luj} = {RviRw} — R[v,w] (|AKN|, 

[^), where the bracket [, ] on C°°{Q*) is induced by the usual commutator (this 
bracket is cahed Lie-Poisson or Berezin-Kirillov-Konstant-Souriau). 

So finding one left-invariant function I2 = Ly commuting with 2H = ■Ci +C2 
will suffice, just take any I3 = i?^. We take the Casimir function F for v, i.e. 
such a function that [-jF] = 0. This function always exists locally since Q has 
odd dimension, but globally it can have singularities. 

In the case of so(3) we have smooth Casimir function F ~ ef + + e^. For 
sl2(K) F — 4eie2 — eg. For the Heisenberg group F — 63. For the solvable 
group we consider [-,63] as a linear vector field on the plane E^(ei, 62). We take 
F to be an integral of this vector field. For example for the case (3a) we take 

Note that these functions F have singularities along some axes in the solvable 
case. To overcome this difficulty we use the lemma. 

Lemma 6. If 3- dim Lie group G is not semisimple, then for every left- invariant 
H there exist two commuting right-invariant integrals I2,l3- 

Actually due to Bianchi classification we can always find a two-dimensional 
commutative subalgebra {v, w) <Z Q in non-semisimple case ((ei, 62) for solvable 
algebras). Now it's easy to check that in every case of non-holonomic metric we 
obtain 3 functionally independent a.e. integrals. □ 

Note on the non-holonomic flovifs on SL2(M). The two contact struc- 



tures on SL2(Z) have nice geometric interpretations ([Kr|). It is well known 
that SL2(Z)/{±1} ~ ST*L^ for the spherical bundle of the Lobachevskii plane 
and since every surface Af^ of genus 5 > 1 can be obtained as quotient of 
by a discrete subgroup we end up with two contact structures on ST*M'^ of 
which one is the standard and the other is the connection form associated with 
a metric of constant negative curvature. 

Note that the metric descends to ST*M'^ and for the second structure even to 
Af^. However we lose integrals: On ST*L'^ we have all two additional integrals, 
on ST*M'^ we loose one right-invariant and under projection to NP we loose 



the other integral. This is consistent with the Kozlov theorem [Ko| according to 
which the geodesic flow on NP with <? > 1 handles is not analytically integrable. 
Still we see there is an integrable lift of the flow. 



5 Non-holonomic metrics on suspensions Ma 

(A) Construction of suh-Riemannian structure. 

Let Ma = T'^xaS^ be quotient of the cylinder C = T'^{ipi,ip2) x R^{ip3) 
w.r.t. the free action given by the map A((pi, (p2, ^3) — {A{ipi, 1P2), ^3 + 1). 

1°) Consider at first the case of semisimple A E SL2(Z) with real eigenvalues 
A > 1,A~^ G (0,1)- Let r]i,ri2 G T^(T'^) be the eigenvectors (the case A = 1 
was considered in example and the case A < will follow similarly because 
the distribution 11 has no canonical orientation and so the change rjj 1-^ —rjj 
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preserves the metrics below). The action on the basis is the following A : 

We define the contact structure on C by the vector fields ^1 = e~^'^^"^''rii + 
ginA <^3j^^^ ^2 = d^p^: n = (^1,^2)- Moreover we fix non-holonomic metric g by 
requiring that ^1,^2 is an orthonormal basis. 

Lemma 7. 11 and g are invariant w.r.t the action A. □ 



Lemma 8 . The distribution H is contact. 
Proof. 

[6,6] =6 = -lnAe-'"^-'^='»7i +lnAe'"^"^^r?2. 
But 6 ^ n because the matrix 



g-lnA-<^3 glnA-¥'3 

-lnAe-^°^-^=* InAe'"^-^!* 



is Vandermonde and hence nondegenerate for In A 7^ 0. □ 
Summarizing we get a sub-Riemannian structure (11,5) on Ma- 

2°) Now let the eigenvalues of A be A, A = 6^^*=. Since the matrix is integer- 
27r 

valued we have 6k = = 3,4, 6 (modulo the cases we have considered). 

Thus A has the matrix ( sm9k\ ^ ^^^.^ 771,772 of TJT^). We 

\^-sm6'fe cos6'fey ' 

define ^1 = cos{dk(P3)vi - sm{ek(p3)il2, 6 = ^^3 and as before fl = (6,6), 
9 — (C*)^ + ($2)^- easily checks that (fl, 5) is invariant, non-holonomic and 
hence defines a sub-Riemannian structure on Ma- 

3°) Consider finally the Jordan box, i.e. A is conjugated to |^ . Let 

771,772 be the corresponding basis, A : (771,7/2) ^-* (771,771+772). We set ^1 = 
cos(27r<^3)77i + sin(27r(/53)(772 — </'377i), 6 = (^<P3 and proceed as before to get 
sub-Riemannian structure (H, g) on M^. 



(B) C^-integrability. 

The Ifamiltonian of the sub-Riemannian flow is if = 5K1 +^2]- Let p',p",P3 
be momenta dual to vectors r]i,r]2,d^^ from C (i.e. evaluation matrix 77(77) is 

£^3x3)- 

This Hamiltonian is A-invariant, so that it descends to Ma- Since H does 
not depend on (pi,ip2 the functions pi,P2 and hence p',p" are integrals on C. 
However these functions are not yl-invariant. We get invariants from their com- 
bination. 

1°) In this case we can write 
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Here the functions transform by the rule A : {p',p",P3) (A ^p' , \p" ,p3). 
Thus we have the following integrals: Ii = H and 

h^p'p", /3^sin(2.i^)e-(^V')-, (7) 

which are invariant and functionally independent a.e. Thus we conclude C°°- 
integrability of sub-Riemannian structure (n,g) on Ma- 

2°) In this case the action is the rotation by 9k and so the additional integrals 

are: 

i2 = {pr + {p'T, h = Mp' + ip'T- (8) 

Note that in this case the integrals are analytic. 

3°) The action is A : {p',p",P3) >^ {p' — p" ,p" ,p3). Therefore the additional 
integrals are: 




This is again the case of C°°, not C"^-integrals. 

(C) Obstructions to "nice" integrability. 
Here we generalize Taimanov's result. 

Definition. We say sub-Riemannian structure (Af'",H, 17) is geometrically 
simple if Q — {H = 1/2} contains a closed nowhere dense invariant subset 
r such that: 

1. Q\T ^ Uj=i J where Uj are open and path-connected. 

2. Each Uj is foliated by Liouville tori, Uj ~ x Z)™"!. 

3'. Let Q be the one-point compactification along fibers of the projection p : 
Q M . Then for every x Cz Q there exists an arbitrary small neighbor- 
hood Wx such that Wx \ T has a finite number of path- connected compo- 
nents. 

The condition 3') is a version of the corresponding 3) from for non- 
compact Q. In fact both conditions are equivalent to the following: 

3". For every x ^ M there exist arbitrary small neighborhoods Ux CI Vx such 
that\7To{[p-HVx\Ux)] u [p~\Ux)\r])\<^. 

We remark also that the set F is usually bigger than the bifurcational set S C Q 
in Liouville integrable case. 

Theorem 9. Let a sub-Riemannian structure {Tl,g) on a closed manifold 
be Liouville integrable and geometrically simple. Then: 
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• The group 7ri(M) is almost commutative, i.e. it contains a commutative 
subgroup of finite index. 

• The first Betti number dimHi{M"^;R) < m. 

Proof. Actually one take a graph with vertices numerated by the com- 
ponents of p~^{Vx \ Ux) U {p^^{Ux) \ r). This number is finite due to 3"). 
Applying sub-Riemannian Hopf-Rinow theorem pR| ] we get labeled directed 
edges of this graph with labels s E Tri{M,x) realized as geodesic loops lifted to 
p~^(yx) n (Q \ r). The rest of Taimanov's proof Q goes now unchanged. □ 

We call the set of a.e. independent involutive integrals Ii — H, I2, . ■ . , In 
of Hamiltonian system on VF^" geometrically simple on Qc = {H = c} if the 
bifurcation sets S are diffeomorphic on Qc' for close c' w c and the stratification 
of Qc by singularities of the moment map Qc ^ x 1-^ {l2{x) . . . , In{x)) can be 
extended to stratification Qc — US" with F = Uj>oS" giving the conditions 
1-3) of the previous definition. Certainly if sub-Riemannian fiow possesses a 
geometrically simple involutive set of integrals, it is geometrically simple. 

For compact Q (Riemannian case) analytic involutive set of integrals is geo- 
metrically simple Q. In sub-Riemannian case non-compactness makes difficul- 
ties for application of the standard theorems about sub-analytic (constructive) 
sets. Yet known integrable cases are geometrically simple. Actually if the inte- 
grals are C"^ and polynomial by momenta (in this case we have Liouville tori), 
the set is geometrically simple. This is because singularities of (/,) give 

stratification of Q at infinity, while the finite (compact) part can be additionally 
substratified (and even made into simplical decomposition as in |t|). 

Corollary 2. Since t:i{AIa) is not almost co mmu tative for A with real eigenval- 



ues different from ±1 and for the Jordan box , GK l, the sub-Riemannian geodesic 



flow is not Liouville integrable with geometrically simple set of integrals. □ 

(D) Calculation of topological entropy. 

Calculations of topological entropy in the case of rotation and Jordan box 
are easy. So we consider only semi-simple case. 

First note that the following submanifold {p' = 0,p" = 0,p3 = 1} of T*C 
is invariant under the transformation A. So this 3-dimensional submanifold 
descends to T*Ma and belongs to the hypersurface {H = 1/2}. Moreover it is 
invariant under the geodesic fiow and the fiow moves with the unit speed along 
9i^3 . So fixing ip3 = we get the torus and the Poincare return map of the 
geodesic flow induces the "Anosov map" T^ with matrix A. Since the 

entropy of a system is not less than entropy of any subsystem, we conclude: 
/itop(sgradiI)>ln|Sp+(A)|. 

Now let us prove the inverse. For this we study bifurcations of (/i ,12,13). Let 
us note that the integrals (M) coincide with the integrals of Bolsinov-Taimanov 



[BT2[, while the Hamiltonians are different. The bifurcation set in the Rieman- 
nian case is twice bigger than in the sub-Riemannian one, but still the systems 
behave similarly. 
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It is easy to see that rk{ci/i , d/2 , dis} = 3 along a trajectory if (p')^ + (p")^ 7^ 
and either 7)3 ^ or e'"''' '''''^' ^ g-in-^' Vay a,t least at one (and then any) 
point of it. Then by the Liouville theorem this trajectory lies on an invariant 
torus and the motion is quasi-periodic. 

Let p3 = ^ ^3 = const and p" = e^i^^ -^^y^ {p')^ + {p"f ^ 0. Then 
the point moves quasi-periodically along the fiber-torus C Ma and hence 
again contributes nothing to the topological entropy. 

Now on th e inv ariant set N — {p'p" = 0} our Hamiltonian (^) coincides with 
the one from | BT2| . Therefore the description of the trajectories is the same. In 



particular N ^ N' N" , each of the summands N' = {p' = 0}, N" = {p" = 0} 
being invariant and diffeomorphic to Ma x . The intersection N' n N" has 
two components V ,V" ~ Ma given by {^3 = ±1} such that V is a stable 
manifold for N' and unstable for N" , V" is stable for N" and unstable for N' . 

Moreover the Hamilton equation of the geodesic flow implies that Lyapunov 
exponents on V are exactly ± In A, 0. So by Katok theorem the only invariant 
Borel measure on N can be supported on V' U V" . But the flow on V and V" 
is Anosov with Poincare time-one map A -.T^ . Therefore the variational 

principle yields hto-p — hi A. 

(E) Reeh vector field. 

We prove in this section theorem |[ 

The unique (up to sign) contact form a for sub-Ricmannian structure is given 
by the equality (ia(^i,^2) = —1 or a([Ci,C2]) = 1, where ^1,^2 is an arbitrary 
(7-orthogonal frame from C°° (11) . 

i°) Consider at first the semisimple case ^1 = ''^'^■'^^?7i -I- e'"'*' '''3ry2, ^2 = 
d^^. The Reeb field satisfies Va — w mod 11 for v = [^1,^2] = In A(e^'"'^ '^3?7i — 

glnA-V3^2). 

Lemma 10 . w is a symmetry of 11. 

Proof. Actually = 0, [t;,6] = In^ Ki- □ 

Corollary Z . Va = v is the Reeh vector field. 

Proof. Actually a{v) = 1. □ 
Now we see that the flow of is quasiperiodic on the tori- fibers of Ma. 
Therefore /itop(j^a) vanish. 

2°) The case of rotation is absolutely similar. 

3°) For the Jordan box the field v — [^1,^2] is not a symmetry of 11 and 
hence is not the Reeb field. Calculation shows that the contact form is a = 
(27r-|-sin^(27r(/33)) (sin(27ri^3)77j[' -|- ((^3 sin(27r(/33) - cos(27r(^3)?72)) , where 771, 772 
is the co-basis of T*T^ dual to r?i,772- Now since da = dip^ A (a), the Reeb 
field is Vq, = V + b{ip^)^i and again hto-piva) = 0. 

Remark 4. Though htop{g) = htonig) in theorem it is not a simple combi- 
nation of corollary [J and theorem formula ^) is not working here. 
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6 Reeb-symmetric SR-structures 

Let n be a contact structure on M^"+^ and Va be the Reeb vector field as 
in introduction. 

Definition. Sub-Riemannian structure (H, g) is called Reeb-symmetric if Va 
is a symmetry of g on H: L^^g = 0. In other words Va is a Killing vector field. 

For instance the Heisenberg SR-structure (example ^) is of this kind. 
Let us define the function Ig{p) = (p, Va)"^ (or just i/^) on T*M . We will call it 
Reeh momentum. 

Proposition 11. (n,g) is Reeh symmetric iff the momentum Ig is the first 
integral of the sub-Riemannian geodesic flow. 

Proof. Actually locally in Darboux coordinates: a — dx ~ X^iLi ^i'^Vi ^^'^ 
Va = dx- Now it's obvious that L^^g = and {Ig, H} — are equivalent to the 
condition that g does not depend on x [lifted from some go on ]R^"(2/, z)]. □ 

Let us construct by g and i^a the Ricmannian metric g on M as in introduc- 
tion. Then the Ricmannian and sub-Riemannian Hamiltonians are connected 
by the formula 

Hg=Hg+Ig. (9) 

Corollary 4 . Reeb momentum Ig is an integral of Ricmannian geodesic flow if 
and only if it is an integral for sub-Riemannian one. □ 

Therefore one can integrate instead of g, the metric g, finding involutive set 
of integrals starting from Ji = Hg and I2 = Ig- 

Remark 5. Let us note that the metrics which equal g onll and l/\/t on 
orthogonal Va, satisfy: lim^distg(t) = distg - sub-Riemannian distance. The 

corresponding Hamiltonians Hgit) — Hg + tig. 

Let us note also about topological entropy: 

Proposition 12. // {Tl,g) is smooth and Reeb-. symmetric, then vanishing of 
two of three htopisgradHg), htopisgradHg), htop{i^a) implies vanishing of the 
third. Moreover if htopii'a) — 0, then htopig) — htop{g). 

Proof. Actually this follows from (^ and commutativity. □ 

Remark 6 . Note that even though the Reeb field is divergence-free and 
shares some good properties (no attractors/repellers) , we can have htopiva) > 0. 
Examples are provided by geodesic flows on negatively curved manifolds. 
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7 Concluding remark: norm via entropy 



Any integrable Hamiltonian system on (M , w) with integrals (/i, . . . , /„) 
determines ]R"-action with generators <~p^ being shifts along sgrad^(/j). Then 
the formula (^) with h = /itop gives a pseudonorm p = ptop on this R" (of course 
for h = hfj, with Liouville invariant measure /i we have = 0) . 

For instance if the action reduces to the action of the torus T" then ptop = 0. 
This is obvious from the definition (also is seen from Atiyah convexity theorem) . 



Examples of this and | BT2 1 papers provide M^-action with ptop ^ 0. However 



Ptop is not a norm, because /itop(sgrad/) — for any / = A2/2 + A3/3. 

Theorem 13. There is an integrable Hamiltonian system with nondegenerate 
Ptop (norm). 

To explain the result consider a piece of bifurcation diagram for the mo- 
mentum map (/i,/2,/3) in the semisimple case that is pictured in Fig. 1. The 
central vertical line corresponds to the most complicated singularity {ps = ±1} 
and in the preimage of each of its point we have 2 copies of M^. Two curved 
boundary surfaces correspond to bifurcation diagram either, while 4 other plane 
pieces of the boundary consist of regular points and appear due to the size of a 
neighborhood. 



Fig. 1 

Now let's take 3 mutually non-intersecting such pieces in the afRne such 
that their central axes are pair-wise non-parallel (we can direct them along axes 
hjhth)- We connect their regular points by small tubes as shown in Fig. 2. 
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Fig. 2 

This is equivalent to connecting preimages of pieces by the tubes T^xlx -D^(e). 
The tubes in M.^ carry the affine structure, that permits by Arnold-LiouviUe 
theorem to define symplectic form and a Hamiltonian on the obtained manifold 
which is integrable with the Hamiltonian action of given by {Ii,l2, h) and 
such that ptop is a norm on this R^. 
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